Boson-fermion mixtures inside an elongated 
cigar-shaped trap 

Z. Akdeniz 1 ' 2 , P. Vignolo 1 ^: and M. P. Tosi 1 

1 NEST-INFM and Classe di Scicnze, Scuola Normale Superiore, 1-56126, Pisa, Italy 
2 Department of Physics, Istanbul University, Istanbul, Turkey 

Abstract. We present mean- field calculations of the equilibrium state in a gaseous 
mixture of bosonic and spin-polarized fermionic atoms with repulsive or attractive 
interspecies interactions, confined inside a cigar-shaped trap under conditions such 
that the radial thickness of the two atomic clouds is approaching the magnitude of 
the s-wave scattering lengths. In this regime the kinetic pressure of the fermionic 
component is dominant. Full demixing under repulsive boson-fermion interactions can 
occur only when the number of fermions in the trap is below a threshold, and collapse 
under attractive interactions is suppressed within the range of validity of the mean- 
field model. Specific numerical illustrations are given for values of system parameters 
obtaining in 7 Li- 6 Li clouds. 
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1. Introduction 

Dilute boson-fermion mixtures have been studied in several experiments by trapping and 
cooling gases of mixed alkali-atom isotopes [H 121 El IH El E] • The boson-fermion coupling 
strongly affects the equilibrium properties of the mixture and can lead to quantum 
phase transitions (for a review see [Zj). In particular, boson-fermion repulsions in three- 
dimensional (3D) clouds can induce spatial demixing of the two components when the 
interaction energy overcomes the kinetic and confinement energies 8J. This occurs in 3D 
if the fermion density is high or if the boson-fermion repulsion is strong, as shown in a 
phase diagram obtained by imposing pressure balance and equality of chemical potentials 
in coexisting phases jUJ EI] • Although spatial demixing has not yet been experimentally 
observed, the experiments of Schreck et al. jTj on 7 Li- 6 Li mixtures inside an elongated 
3D trap appear to be close to the onset of a demixed state [TT]. Collapse driven by 
increasing the number of bosons in a 3D trap has instead been observed in 87 Rb- 40 K 
mixtures [Sj, where the interspecies interactions are attractive. 

On approaching a quasi-two-dimensional (Q2D) situation, the Fermi energy starts 
to increase linearly with the number of fermions per unit area, so that the linear stability 
condition of the mixture becomes independent of the fermion density and involves only 
the scattering lengths and the transverse width of the cloud. In this regime spatial 
demixing of a 7 Li- 6 Li gas is predicted to occur by simply squeezing down the axial 
thickness of the atomic clouds [T2j . 

The analysis of Viverit et al. [H] for homogeneos boson-fermion matter has been 
extended by Das to a one- dimensional (ID) homogeneous mixture. In the ID case 
the Fermi energy is proportional to the square of the number of fermions per unit 
length, so that the mixed state becomes unstable at low linear fermion densities. The 
dependence of the kinetic pressure of the Fermi gas on the fermion density becomes 
stronger with decreasing dimensionality, and in ID a number of fermions above a 
threshold will want to occupy the whole available space irrespectively of the repulsions 
exerted by the bosons. Similarly, it is easily shown that a ID homogeneous mixture 
with attractive boson-fermion interactions is stabilized against collapse by the Fermi 
pressure if the linear density of fermions is above the same threshold. 

In this paper we examine the instabilities of a dilute boson-fermion mixture trapped 
inside a strongly elongated harmonic confinement. We give the general conditions under 
which various forms of demixing can occur in a Q1D confined geometry in terms of 
the physical parameters of the system and draw a phase diagram at zero temperature 
in a plane defined by two scaling parameters. We support our analysis by extensive 
mean-field calculations of the equilibrium density profiles in a 7 Li- 6 Li mixture. We find 
that in the case of interspecies repulsions the stability condition which applies in the 
macroscopic limit is still approximately valid for full demixing in harmonic confinement. 
We next turn to the case of attractive interspecies interactions, where we search for a 
collapse instability by a variational approach patterned after the work carried out by 
Miyakawa et al. [TJ] on a 3D mixture. We find that our mean-field model does not allow 
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collapse in a confined Q1D geometry, but only predicts that the density profiles keep 
narrowing with increasing attractive coupling strength until the validity of the model 
is lost. These predictions of the variational analysis are again supported by numerical 
illustrations of density profiles. 

The paper is organized as follows. In Sec. El we describe the mean-field model 
and report the stability condition for a macroscopic mixture in ID and Q1D geometry. 
Section El discusses the conditions under which various forms of demixing occur in a Q1D 
harmonic confinement and shows some illustrative configurations of the density profiles 
in the gaseous cloud in the phase-separated regime. The variational analysis and some 
illustrative density profiles for the case of interspecies attractions are given in Sec. |U 
Finally, Sec. El gives a summary and some concluding remarks. 



2. The model and the phase diagram 

The atomic clouds are trapped in cigar-shaped potentials given by 

V b e f(x,y,z) = m bJ ul {bJ) {x 2 + y 2 )/2 + m bJ uj 2 zib j ) z 2 /2 = U bJ (x,y) + V bJ (z) (1) 

where m b j are the atomic masses and uJ z (b,f) UJ ±(bj) the trap frequencies. We assume 
repulsive interactions between the bosons and focus on the case where the trap is 
elongated enough that the dimensions of the atomic clouds in the radial plane, which 
are of the order of the radial harmonic-oscillator lengths a±(b,f) = (h/ m bj UJ ±(bj)) 1 ^ 2 , 
are comparable to the magnitude of the 3D boson-boson and boson-fermion scattering 
lengths, a bb and \a b f\. In this regime we can study the equilibrium properties of the 
mixture at essentially zero temperature (T ~ 0.02 Tp) in terms of the particle density 
profiles along the z axis, which are n c (z) for the Bose- Einstein condensate and nj(z) for 
the fermions. The profiles are evaluated using the Thomas-Fermi approximation for the 
condensate and the Hartree-Fock approximation for the spin-polarized fermions. Here 
and in the following we shall take a± b = a±f (= a±, say). 

The Thomas-Fermi approximation assumes that the number of condensed bosons 
is large enough that the kinetic energy term in the Gross-Pitaevskii equation may be 
neglected [T3]. It yields 

n c (z) = \ji b - V b (z) - g bf n f (z)]/g bb (2) 

for positive values of the function in the square brackets, and zero otherwise. Here, 
fi b is the chemical potential of the bosons. This mean-field model is valid when the 
high-diluteness condition n c a bb -C 1 holds and the temperature is outside the critical 
region. If the conditions a± > (a bb , \a b f\) are fulfilled, the atoms experience collisions in 
three dimensions and the coupling constants in Eq. can be written in terms of those 
for a 3D cloud as [THl HZj 

_ _ 9bb _ 9 b f r „v 

with g bb = 47ih 2 a bb /m b , g b j? = 27ih 2 a b f jm r and m r = m b mf/(m b + uif). 
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n f (z) 



(4) 



The Hartree-Fock approximation [T%1 IT^l I2U] treats the fermion cloud as an ideal 
gas subject to an effective external potential, that is 

/^M(4 + ^'w-«)H +1 P 

where ///is the chemical potential of the fermions and 

V eff (z) = V f (z)+g bf n c (z). (5) 

The fermionic component has been taken as a dilute spin-polarized gas, for which the 
fermion-fermion s-wave scattering processes are inhibited by the Pauli principle and p- 
wave scattering is negligibly small [21] ■ In the mixed regime the Fermi wave number in 
the axial direction and the condensate density should be smaller than l/|a&/|, but this 
is not a constraint in the regime of demixing where the boson-fermion overlap is rapidly 
dropping. On the contrary, for large attractive interactions this condition, together 
with the diluteness condition n c a^ <C 1, may not be fulfilled and this sets the range of 
validity of our results. 

The chemical potentials fi^j characterise the system in the grand-canonical 
ensemble and are determined by requiring that the integrals of the densities over the z 
axis should be equal to the average numbers N b and Nf of particles. The presence of 
a bosonic thermal cloud can be treated by a similar Hartree-Fock approximation [11], 
but it has quite negligible effects at the temperatures of present interest. 

As the boson-fermion coupling increases, the mixture can become unstable against 
demixing (in the case g b f > 0) or collapse (in the case g b f < 0). In the macroscopic ID 
limit at given particle densities a linear stability analysis based on a mean-field energy 
functional predicts that the locations for demixing and collapse coincide and that the 
gaseous cloud is stable if the usual condition 

9bb9ff -9b f >0 (6) 

is fulfilled 0, with gjf = h 2 7T 2 nf/mf playing the role of an effective fermion-fermion 
repulsion due to the Pauli pressure in the Fermi gas. Since #// increases with rif, the 
mixed state is stable at given boson-fermion attractive or repulsive coupling strength if 
the linear density of fermions is above a threshold 13J. 

The stability condition in Eq. © can be extended to the Q1D geometry for large 
numbers of particles and if the boson and fermion densities vary smoothly (for the sake 
of brevity in the following we will refer to this regime as the "Q1D macroscopic limit"). 
By using Eq.© for the boson-boson and boson-fermion couplings and by taking the 
value of rif as the linear fermion density at the centre of the trap, Eq. @ is rewritten 
as 

\at\<(*mi k faa ) 1/2 (7 ) 

a± \mj,mf ) 

where kj = (2Nf) 1 ^ 2 /a z f is the Fermi wave number under harmonic confinement at the 
centre of the trap, with a z f the axial oscillator length for the fermions. In Fig. E]we 
have plotted Eq. (J7|) in the plane kfQbb} for the case of a 7 Li- 6 Li mixture. 
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Figure 1. Phase diagram of a Q1D 7 Li- 6 Li mixture at T = 0, in the plane defined 
by the dimensionless parameters atbkf and \at,f\/a±. The symbols at \abf\/a± = 0.5 
correspond to the parameters of the configurations shown in Figs. and EJ (□) 
N f = 3 x 10 5 , (O) N f = 8,5 x 10 4 , (+) N f = 2 x 10 4 , (x) N f = 3740. 

The effects of finite size and dishomogeneity of the gaseous clouds are not taken 
into account in Eq. (JZj), with the consequence that the resulting phase diagram in Fig. 
[0 shows a sharp transition as in the macroscopic limit. Transitions in a mesoscopic 
cloud under confinement are instead spread out and several alternative definitions of 
their location can therefore be given. These are discussed in the next two sections. We 
shall for the sake of brevity continue to use expressions such as "phase transitions" and 
"phase diagram" when needed. 

3. Spatial demixing 

In the calculations that we report in the following we have used values of the trap 
frequencies appropriate to the experiments in Paris on the 7 Li- 6 Li mixture pQ, that is 
u z b/2n = 4000 Hz and uj z f/2-n = 3520 Hz. We also take a&& = 5.1 ao, a being the 
Bohr radius, as appropriate to 7 Li- 7 Li s-wave scattering. We first discuss in this Section 
three alternative locations of the demixing point in a mesoscopic trapped cloud, that 
we denote as partial, dynamical, and full demixing. Their definition and the derivation 
of simple analytical expressions in the Q1D model are given below. 

3.1. Partial demixing 

The interaction energy E int between the boson and fermion clouds initially grows as the 
boson-fermion coupling is increased, but in the demixing regime reaches a maximum and 
then falls off as the overlap between the two separating clouds diminishes. We locate 
partial demixing at the maximum of Ei nt , that we calculate from the density profiles 
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E in t = 9bf / dzn c (z)n f (z). 



The behaviour of the interaction energy as a function of a&/ at a± 
various numbers of bosons and fermions, is shown in Fig. |21 
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Figure 2. Boson-fermion interaction energy E int (in units of Nhoj z i,, with N = 
Nb + Nf), as a function of a b f (in units of aj_) for a bb = 5.1 ciq. Left panel: 
N b = N f = 10 4 . Right panel: N b = 100 and N f = 19900. The arrows indicate 
the locations of demixing estimated from Eqs. 10, (|13|) . and (|14f) . 



An approximate analytic formula for the location of partial demixing as a function 
of the system parameters in the Q1D model can be obtained from the condition 
dEi n t/dgbf = which leads to the relation n c (z)rif(z) — nj(z)/gbb — n l( z )/9ff = 0. By 
using the estimate n c j rs N^j / (2Rbj) with the values of the cloud radii in the absence 
of boson-fermion interactions, Rf = (2Nf) 1 ^ 2 a z f and Rb = [SNbgbbdzb/ 1 {2fvjj z b)} l ^a z b, 
where a z (b,f) — (^-/ m 6,/^(fe,/)) 1//2 ; the location of the maximum in the boson-fermion 
interaction energy as a function of a,bf/a± is found to lie at 

-l 



a bf 

a -L part 

where we have defined 



Ci kfd bb + 



C-2 



k 2 n 2 



ci 



C-2 



m b a zb ( 3a bb a 



zb 



AN, 



2/3 



.a± \ 



AT 2 / 3 

N h m f abb 
it 2 m r a z b V 3a±a z b 



a bb 



1/3 



1/3 



(9) 



(10) 



The prediction obtained from Eq. Q is shown in Fig. EJby solid arrows. For the case 
Nb = Nf = 10 4 , there clearly is reasonable agreement between the analytical estimate 
and the numerical results. For a low number of bosons and a large number of fermions 
the condition given in Eq. Q is dominated by the term involving c\. In this limit, the 
bosonic density is largely deformed by the pressure of the fermions and the prediction 
obtained from Eq. Q overestimates the position of the maximum of the interaction 
energy (see right panel in Fig. EJ). 
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3.2. Dynamical demixing 



H (12) 



At some value of the boson-fermion coupling the fermion density vanishes at the centre 
of the trap. This occurs when 

9bf_ _ __ 

9bb dyn 

We denote this point as the dynamical location of demixing in a mesoscopic cloud, since 
we expect a sharp upturn of the low-lying fermion-like collective mode frequencies to 
occur at this point, where the topology of the cloud changes. This has been shown to the 
case for both collisional and collisionless excitations in a mixture under 3D confinement 



If we insert the chemical potentials for ideal-gas clouds in Eq. (112)1 . we obtain an 
approximate expression for the location of dynamical demixing, 

/ 2 \ 1/2 

k f a bb . 13 

\m b mf J 



a bf 
a± 



dyn 

The prediction obtained from Eq. (|13|) is indicated in Fig. El by dot-dashed arrows. 
For both small and large values of the number of bosons, the dynamical condition 
for demixing is in good agreement with the numerical results from the boson-fermion 
interaction energy. 

3.3. Full demixing 

The point of full demixing is reached when the boson-fermion overlap becomes negligible 
as in a phase transition occuring in a macroscopic cloud. If we use the stability criterion 
given in Eq. ((Zj) for the macroscopic Q1D case, the condition for full phase separation 
at T = inside a strongly elongated trap is given by 

2 \ V2 

r k f a bh (14) 



a bf 



( 27r?;. ! 



full \m h m f 

or a b f\f u u ~ \f 7 2,a b f\ ( iyn- The prediction obtained from Eq. (|THl is indicated in Fig. El 
by short-dashed arrows. Analogously to the partial demixing condition there is better 
agreement with the numerical results on the boson-fermion interaction energy when the 
number of bosons is relative large. 

In summary, in Q1D a small number of fermions allows demixing, which is in 
contrast to the case for full phase separation in 3D where the critical value of the boson- 
fermion scattering length decreases by increasing the number of fermions [5J , and to 
the Q2D case where the transition is controlled only by the geometrical parameters of 
the trap and not by the number of atoms [T2j. 

The critical coupling strength for both dynamical and full demixing under harmonic 
confinement moves upwards with the law g b f oc N l ^ A as Nf increases, till the mixture 
becomes stabilized by the kinetic pressure of the fermions. 
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Figure 3. Density profiles of the condensate (left panel) and of the fermions (right 
panel) at abf/a± = 0.5 for various numbers of fermion atoms: Nf = 3 x 10 5 (solid 
line), N f = 8.5 x 10 4 (dashed line), Nf = 2x 10 4 (short-dashed line), and Nf = 3740 
(dotted-dashed line). Notice the difference of the horizontal scales in the two panels. 



3.4- Density profiles at full demixing 

With the aim of verifying the validity of the condition of full component separation in 
Eq. ()14J) , we have evaluated the boson and fermion density profiles by varying the boson- 
fermion scattering length as might be attained experimentally by exploiting optically or 
magnetically induced Feshbach resonances [23]) or by varying Nf i.e. kf at fixed radial 
size a_i_ = 38 do- As illustrative examples we show in Fig. El the density profiles of the 
lowest-energy configurations indicated by the symbols in Fig. ^ In agreement with the 
location of the configurations in the phase diagram, density profiles with a small number 
of fermions are spatially demixed, whereas stability is reached on increasing the number 
of fermions. 

In the region of instability we have also found in addition to the thermodynamically 
stable configurations, shown in Fig. El various metastable "exotic" configurations for the 
demixed clouds. Some illustrative examples of such metastable density profiles, which 
break axial symmetry and/or consist of alternating slices of bosons and fermions, are 
shown in Fig. 0] 



4. Collapse 



In the case of attractive boson-fermion interactions we have studied the instabilities of 
the Q1D harmonically trapped mixture by a variational method. Following Miyakawa et 
al. [14 1 . wc write the total energy of the g clS clS db function of the boson order parameter 
$(2), with |$(z)| 2 = n c (z), and of the fermion density distribution n/(z) as 

h 2 



E[9{z),n f (z)] 



dz 



2mb 



w (z ) 2 + v b (z m (z) 



+ f |$(z)| 4 



dz 



6m 



n)(z) + V f (z)n f (z) + g bf n f {z)\^{z)[ 



(15) 
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Figure 4. Metastable density profiles of fermions (solid line) and bosons (dashed 
line) at full demixing, for a Q1D mixture of Nb = 10 4 bosons and Nf = 4 x 10 4 
fermions with a b b = 5.1 do, a b f = 19ao, and a b f/a± = 0.5. The various configurations 
correspond to an excess energy AE/E over the thermodynamically stable demixed 
configuration given by AE/E ~ 0.2% (left-top panel), AE/E ~ 1% (right-top panel), 
AE/E ~ 1.6% (left-bottom panel), and AE/E ~ 5% (right-bottom panel). 



We neglect the kinetic energy of the bosons and make a Thomas-Fermi Ansatz for the 
boson order parameter, 



~2 \ 1 V 2 

R 2 ) 



(16) 



where z = z/a zb and the variational parameter R is the dimensionless Thomas- Fermi 
radius of the condensate. Minimization of the free energy E[$(z), rif(z)] — fif(R)nf(z) 
yields the fermion density distribution 

_ £2 \ l ' 2 
X 2 ) ' 



n f (z;R) 



2Nf 



iv\a 



zf 



17) 



where A ~ N)' 2 [%R 3 /{Zg^N^] 1 ^ is the t urning point for the fermions at large 
interaction strength g b f = \g b f\/(a zb hu> z f). Finally assuming that A > R, i.e. as long as 
the fermion cloud is much broader than the bosonic one, the total energy E(R) becomes 



E(R) 



zb 



/3m f uj\ 1/2 
V 1m h u\ ) 



(9bfN b ) 



N 2 

1/2 iV / 
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4R 3 / 2 10Ra zb hv b 
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3g bf \^5N^Nfg bf 
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N h 



8 J 4nR 3 / 4 a zf huj b y/24u} b (g bf N b y/ 



_ R V2 + ^L R 2 (18) 
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Figure 5. Energy E(R) (in units of Hu>zb) of a boson-fermion cloud in Q1D harmonic 
confinement, as a function of the variational parameter R for Nb = Nf = 10 4 . The 
curves correspond to a,bf/a± = —0.26 (solid line), —0.5 (dashed line), and —1 (short- 
dashed line). 




Figure 6. Density profiles of the condensate (left panel) and of the fermions (right 
panel) at a,bf/a± = —0.5 for various numbers of fermion atoms: Nf — 3 x 10 5 (solid 
line), Nf — 8.5 x 10 4 (long-dashed line), Nf = 2 x 10 4 (short-dashed line), and 
Nf = 3740 (dotted-dashed line). 



The validity of this expression is limited by the diluteness condition 2Ra zb 3> 
N b max(a bb , \ a bf \). 

Contrary to the 3D case, in Q1D the boson-fermion attraction energy does not 
modify the leading term of E(R) in the small-i? region, so that E(R) remains positive 
and for R — > (see Fig. EJ. The mean- field model only leads to narrower and narrower 
density profiles, until pressure is relieved by the expulsion of atoms from the trap through 
three-body processes allowing collapse to proceed further. Some illustrations of mean- 
field density profiles in correspondence to the symbols shown in Fig. is given in Fig. 
13 The condensate density increases at the centre of the trap with a decrease in the 
number of fermions, which is thus seen to favour collapse. 
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In summary, we have evaluated a mean-field description of a boson-fermion mixture 
confined inside a cigar-shaped trap, such that the scattering events can still be considered 
as three-dimensional but nevertheless affected by the radial confinement. We have 
focussed on the equilibrium properties of the mixture and in the macroscopic limit 
we have given a universal phase diagram for the transitions to demixing or collapse 
expressed in terms of scaling parameters of the system. 

We have studied the boson-fermion interaction energy in a mesoscopic cloud as a 
function of a boson-fermion repulsive coupling and have given approximate analytical 
expressions identifying three critical regimes: (i) partial demixing where the boson- 
fermion interaction energy attains maximum value from a balance between increasing 
interactions and diminishing overlap; (ii) dynamical demixing where the fermionic 
density drops to zero at the centre of the trap and a sharp dynamical signature of 
demixing may be expected; and (iii) full demixing where the boson-fermion overlap is 
negligible as in the macroscopic limit. It is remarkable that in this model full demixing 
can be reached by simply decreasing the number of trapped fermions. In the region of 
phase separation we have found several metastable configurations for the demixed cloud, 
having various topologies but lying at higher energy above the stable configuration which 
is composed of a core of condensed bosons surrounded by fermions. 

When the boson-fermion interaction is instead attractive, a Q1D mixture reaches 
collapse in the macroscopic limit under the same conditions as it attains full demixing 
for repulsive coupling, i.e. upon lowering the fermion density. While a mean-field 
treatment is clearly unable to describe the full path to collapse, we have seen that 
it supports the above view. It would be interesting to study the Q1D mixture with 
strong attractive interactions by using models beyond mean-field and with inclusion of 
three-body collisions. 
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acknowledges support from TUBITAK and from the Research Fund of the University 
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